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ABSTRACT 

We present HerschelSV^RE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar 
dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, sup- 
porting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm 
temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investiga- 
tions. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of 
environmental effects on galaxy evolution. 
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1. Introduction 

Clusters of galaxies are extremely hostile environments for star- 
forming galaxies. A plethora of observations and numerical sim- 
ulations have revealed how gravitational a nd hydrodynamical 
interactions can affect cluster spirals (e.g., iBoselli & Gavazzil 
l2006h : stars and gas can be stripped, galaxy morphologies can 
be changed and star formation activity can be enhanced and/or 
quenched. However, still very little is known about the effects 
of the environment on the dust content of cluster galaxies. Since 
the dust is mixed with the interstellar medium (ISM), the com- 
mon expectation is that the environment should affect the dust 
in a similar fashion as the gas. However this hypothesis has 
still to b e confirmed observationally. The launch of Herschel 
(iPilbratt et al. 2010) has opened a new era in the study of dust 
in galaxies. Thanks to its high spatial resolution and sensitivity 
to all dust components, Herschel should be able to determine 
whether dust is stripped from infalling cluster spirals and dis- 
persed into the intra-cluster medium (ICM). 

One of the most dramatic examples of the effects of the en- 
vironment on nearby galaxies is represented by the disturbed 
galaxy NGC4438, in the Virgo cluster. NGC 4438 is a highly 
Hi-deficient early-type spiral showing stellar tails (Kennev et al. 
1995: Bose lli et al.1 l2005l) and extra-planar gas (IVoUmer et al.l 
2009). Numerical simulations show that only a combination of 
ram-pressure stripping and tidal interactions is able to repro- 
duce the dis turbed morphology and kinematics of NGC4438 
(IVollmer et al. 2005). However, while Combes et al. (1988) and 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



IVollmer et all (l2005h invoked a high- velocity (-800 km s"^) 
gravitational inter action w ith the companion early-type galaxy 
NGC4435, Kenne yet al.l (|2008) have revealed the presence of 
a series of HQf-h[NII] filaments connecting NGC4438 to the gi- 
ant elliptical M86 (see Fig.[T]). These new observations favour a 
more complex scenario in whic h NGC4438 has rece ntly (-100 
Myr ago) interacted with M86 (iKennev e t al."'2008'). Although 
the detailed history of NGC4438 is still unclear, its peculiar 
properties make it an ideal target to test the power of Herschel 
in unveiling the effects of the environment on the dust properties 
of cluster galaxies. 

In this Letter we present Herschel-SVW^. (Griffin et~aD 

I 201Q) observations of NGC4438 obtained as part of the Herschel 
Science Demonstration (SD) Phase. The observations of the 
elliptical galaxy M86 are presented in a companion paper 
(Gomez et al. 20 id ) We a ssume for NGC4438 a distance of 17 
Mpc (Gavazz i et aPl 1999b . corresponding to a linear scale of 82 

2. Observations and data reduction 

The region around NGC4438 was observed by the Herschel- 
SPIRE instrument as part of the SD observations for the 
Herschel Reference Survey (Boselli et al. 2010a). Eight pairs 
of cross-linked scan-map observations were carried out over an 
area of ~12'xl2' with a nominal scan speed of 30'7sec. The 
data have been red uced following th e proce dure described in 
iPohlen et al.1 (|2010') and 'Bendo et alj (l2010al) . However, since 
in this case the median baseline of the whole time line (prior 
to the de- striper) left some residual large-scale gradients (see 
iPohlen et aDl2010l) . we preferred an alternative baseline subtrac- 




Fig. 1. Th e NGC4438/M86 re gion as seen by SPIRE. This mosaic is obtained by combining the 250yum maps centered on NG C4438 
and M86 (iGomez et al .1120 id) . The blue contours show the extended Ha+[NII] filaments discovered by lKenney et al.l (l2QQ8h . Note 
that the intensity excess south of NGC4438 (near the edge of the frame) is just a residual of the baseline subtraction. 



tion on a scan by scan basis. A robust linear fit with outlier re- 
jection was applied to the first and last fifty sample points of the 
time line for each bolometer, thus avoiding the galaxies and the 
extended diff'use emission present across the field. The SPIRE 
astronomical calibrati on methods and accuracy are outlined in 
ISwinyard et al.l ([2010). Since at the time of the data reduction 
the SPIRE pipeline used a preliminary flux calibration, we fol- 
lowed the recommendation of the SPIRE Instrument Control 
Center and multiplied the flux densities by 1.02, 1.05, and 0.94 
at 250, 350, and 500 yum, respectivel}Q. The images have flux 
calibration uncertainties of 15% and the rms are -5.5, 5.8, 6.9 
mJy/beam at 250, 350 and 500 yum, respectively. We note that in 
all three bands the noise is dominated by confusion. The astro- 
metric uncertainty is ~2'' and the full widths at half maximum 
of the SPIRE beams are 18.1'', 25.2'', 36.9" at 250, 350 and 500 
yum, respectively. 



3. Results 

In Fig. [T] we show a m osaic of the 250 yum SPIRE maps 
of NGC4438 and M86 (iGomez et al.l 2010). In addition to 
NGC4438, also the companion SO galaxy NGC4435 is clearly 
detec ted in all the three SPIRE bands. This is not surprising 
since iPanuzzo et al.] f2007) showed that this galaxy has just ex- 
perienced a burst of star formation and it is confirmed by the 
fact that the SPIRE flux density ratios (/(250)//(350) -2.9 and 
/(350)//(500) -2.6) are consistent with th e values observed in 
star- forming galaxies (iBoselli et al. ll2010b'). 

More intriguing is the presence of a few bright (i.e., peak 
/(250) >0.05 Jy/beam) extended sources possibly associated 
with NGC4438. The most remarkable feature is the extended 
emission to the N-NW of NGC4435 (Plume in Fig. [B. This 
plume is approximately 8' long and 3' wide and it is detected 
(at least in part) in all the three SPIRE bands. Alt hough its ori- 
in (i.e .. Galactic or extragalactic) is still debated. ICortese et aP 
2010ah have recently shown that the properties of this plume 



are more consistent with Galactic cirrus than tidal debris at the 
distance of the Virgo cluster. 

Two additional elongated structures (El and E2 in Fig. O 
are detected within the optical radius of NGC4438. El appears 
as a tail extending from the bulge of NGC4438. However a care- 
ful comparison of the SPIRE data with Spitzer images reveals 
that it is just the result of blending of point sources visible at 
3.6 and 24 yum, possibly background galaxies. This seems also 
supported by the lack of any counterpart in UV, Hi or HcK-h[NII]. 
Less clear is the origin of E2, which is composed of three dif- 
ferent knots. The south ern and cen t ral kn ots coincide with an 
Hi cloud discovered by iHota et all (12007 ) and the whole fea- 
ture apparentl y follows one of the HQf-F[NII] filaments pointing 
towards M86 (Kenney et al.l 120081) . suggesting that it might be 
associated with NGC4438 (see Fig. O. If so, this would repre- 
sent a unique example of stripped intra-cluster dust. However, 
caution is required before interpreting E2 as a dust stream in 
Virgo. Firstly, the southern knot is off'set (~15" north) from the 
HQf-h[NII] stream, thus it is not clear whether the two features are 
related. Secondly, Spitzer 3.6 and 24 yum point sources are found 
within ~6-9" from each of the three knots, thus we cannot com- 
pletely exclude the possibility that E2 is just the result of blend- 
ing of background sources, like El. So, although E2 is a very 
intriguing system, only future investigations will reveal whether 
it is associated with NGC4438. We note that, in this case, E2 
would represent the only case of submillimetre (submm) emis- 
sion associated with the intra-clu ster Hof-FfNII] streams connect- 
ing M86 and NGC4438 (see also lGomez eTaf . 2010). 

Stronger evidence supporting dust stripping is found in the 
main body of NGC4438. In Fig. [2] we compare the distribution 
of the cold dust as revealed by Herschel with those of stars, 
hot/warm dust and warm ionizecfl cold atomic and molecular 
hydrogen. Interestingly, despite the presence of a young stellar 
population, no submm emission is detected in the northern and 
southern tidal tails. Although the 250-500 yum emission is mainly 
concentrated within the central -4.5 kpc of the galaxy, extended 



^ See also http://herschel.esac. esa.int/SDP-wkshops/presentations7| 
IR/3_Griffin_SPIRE_SDP2009.pdf. 



2 We note that the Hof+LNII] in NGC4438 and the stripped fila- 
ments is likel y not t racing star formation but gas cooling aft er shocks 



ments is likely not t racing star formation but gas cooling ait er 
dKennev et al.lll995l : lMachacek et al.ll2004i : IVollmer et al.ll2009.) . 
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Fig. 2. A panchromatic view of NGC4438. Each column shows the distribution of a different baryonic component. First column: 
stars (GALEX UV, SDSS gri, 2MASS JHK). Second column: cold dust (SPIRE 250, 350, 500 /im)^hird column: warm dust 
(Spitzer 8, 24, 70 jum). Fourth column: cold atomic (Hi, Hota et al. 2007) , molecular (CO. Vollmer et al.ll2005l) and warm ionized 
rHa+rNII1. [Kennev et a02008h hydrogen. The contours show the 250//m emission. Contours levels are 0.03, 0.04, 0.07, 0.1 1, 0.15, 
0.19 Jy beam"^. The features discussed in § 3 are highlighted in the 250 jim map. 



emission is visible on the west side of NGC4438. Particularly 
remarkable are the two bright regions highlighted in Fig. [2l a 
tail (A) extending to the south up to ~9 kpc from the center 
of NGC4438, and an elongated structure (B) ~4.5 kpc from the 
plane of the disk and almost completely detached from the main 
bodyofNGC4438. 

The tail A is detected at 8, 24, 70 //m and in Hi. It extends up 
the southern stellar tail of NGC4438 and is composed of at least 
two disjoint bright knots. Although it follows very closely the 
southernmost HQf+[NII] filament in NGC4438, a careful com- 
parison between the 8, 24 yum and HQf+[NII] images reveals that 
two features do not spatially coincide. The fact that the Hi in cor- 
respondence of tail A has a recessional velocity (-40< V < -h20 
km s"^ , Hota et al. 2 007) significantly low er than the ionized gas 

(V 85 km s"^ iChemin et al.ll2005h suggests that the two 

components might just be projected along the same line-of- sight. 



Although the presence of submm emission up to the stellar tidal 
tail is consistent with a dust stripping scenario, it remains un- 
certain whether tail A is just part of what is left of the disk of 
NGC4438 or a stripped dust tail. 

In comparison, the formation history of knot B appears a lit- 
tle bit clearer. This feature coincides with the extra-planar dust- 
lane visible in optical images a nd it clearly follows the distribu- 
tion of the extra-planar atomic (iHota et al.ii200 7) and molecular 
(Vollmer et al. 2005) hydrogen detected in this region. While the 
extra-planar CO(l-O) emission is mainly segregated around knot 
B, the Hi extends further south following closely the diff'use 250 
yum emission. Although low surface brightness Ha-h[NII] emis- 
sion is observed in correspondence of the diff'use 250 yum to the 
west of NGC4438, only a single Hii region (slightly off'set from 
the submm emission peak) is observed in knot B. Previous inves- 
tigations have shown that the disturbed morphology and kine- 
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matics o f the gas in this region is consistent w i th a strippin g 
scena rio ('Kennev et al.l ll995l:ICombes et al.lll988l: IVollmer et al. 
l2QQ5l) . So, it is likely that we are directly witnessing dust in the 
process of being removed from the disk of NGC4438. Contrary 
to tail A, knot B is not entirely detected by Spitzer. only the sin- 
gle Hii region and low surface brightness emission are visible at 
8 and 24 yum. The integrated /(350)//(24) flux density ratio of 
knot B (~47) is a factor ~4 higher than the value observed in the 
main body of NGC4438, confirming that this feature is missing 
a warm dust component and it is not associated with active star 
formation. This is additionally supported by the reddening of the 
/(250)//(350) ratio (from -2.7 to -2.1) when moving from the 
center of NGC4438 to knot B. 

4. Discussion and conclusions 

The SPIRE data alone are not sufficient to determine whether 
the extra-planar dust is in the process of being removed from 
NGC4438 or, for example, falling back onto the disk. However, 
by combining multiwavelength observations with detailed nu- 
merical simulations, Vollmer et al. (2005, 2009) have shown that 
the distribution and kinematics of the different components of 
the ISM in NGC4438 can only be reproduced via a combination 
of tidal interaction and ram-pressure stripping. Although the de- 
tails of the gravitational interaction are still unclear (i.e., M86 
and/or NGC443 5 „ Vollmer 2009), it appears that strong on-going 
ram-pressure (in addition to tidal forces) is necessary to repro- 
duce the properties of the extra-planar gas component in knot B. 
In this case, whatever the exact mechanisms affecting NGC4438, 
it is likely that at least part of the dust in the west side of the 
galaxy is in the process of being removed from the disk thus 
providing evidence of dust stripping by environmental eff'ects. If 
completely removed, the stripped dust will be dispersed in the 
ICM contributing to its metal enrichment. 

Since in galaxies the dust is associated with the gaseous com- 
ponent of the ISM, it is generally expected that when the gas is 
stripped part of the dust can be removed as well. What is still un- 
clear is how much dust follows the fate of the stripped hydrogen. 
Unfortunately, the SPIRE fluxes alone are not sufficient to accu- 
rately quantify the amount and temperature of the dust in knot 
B. Nevertheless, we can at least try to estimate the dust mass 
by using the ffuxes at 250 and 350 yum, where knot B is clearly 
resolved, and assuming M^ust = (fvD^)/[f<vB(y, T)], where fy is 
the flux density (-0.98 and -0.46 Jy at 250 and 350 yum, respec- 
tively), D is the distance and Ky is the absorption cross section per 
mass of dust. We adopt Ky=4 and 1.9 cm^ g"^ at 250 and 350 yum 
respectively (Draine 2003). The total mass of dust in knot B is in 
the range -2x10^-2x10^ M© for a dust temperature between 10 
and 20 K. In order to determine the total gas mass in the extra- 
planar cloud, we combined the estimates of M(H2)-4.7xlO^ 
Mq (including helium) a nd M(Hi) -1.5x10 ^ M© obtained by 
IVoUmer et aP ([2005) and lHotaetal.1 (l2007h . respectively. The 
total gas-to-dust ratio of knot B is in the range -30-300, not 
significan tly different from the values observed in n earby galax- 
ies (e.g., iDraine et aP l2007l: iGalliano et al.l |2008|) and across 
the disk of st ar- forming spirals (e.g.,'Munoz-Mat eos et al.ll2009l: 
iBendo et al.^ 2010b; Pohlenetal. 2010). This seems to support 
a scenario in which a cloud of gas does not lose a significant 
amount of its heavy elements when removed from the disk. It is 
in fact quite remarkable how well the cold dust follows the spa- 
tial distribution of the cold hydrogen across the whole galaxy, 
even in such a highly perturbed system like NGC4438. 

Interestingly, the high degree of spatial correlation between 
cold dust and cold hydrogen in NGC4438 is different from what 



is observed by Herschel in M86 (iGomez et al .1120101) . In M86 
the dust is mainly associated with the HQf-F[NII] emitting gas 
and only mildly correlated with the cold atomic component. This 
is intriguing and might suggest different dust properties (e.g., 
stripping mechanism, heating source, etc.) in the two galaxies. 

Finally, we note that it is unclear whether the extra-planar 
dust is missing a hot component just because the interstellar ra- 
diation field outside the plane is not strong enough to keep it 
hot, or if the cloud was already dominated by a cold dust com- 
ponent while in the disk. Although the /(250)//(350) flux ratio 
is lower than the value observed in the central part of NGC4438, 
it is not significantly different fr om what is obs erved in the outer 
parts of the disk of M81 (Bendo et al. 11201 Oal) . M99 and MlOO 
(lPohlenetal.ll2010h . 

In summary, in this paper we provide evidence for dust strip- 
ping by environmental effects in NGC4438. The high spatial res- 
olution and sensitivity to cold dust of the SPIRE camera allowed 
us to discover an extra-planar cold dust component completely 
missed by previous far-infrared surveys. The strong spatial cor- 
relation between cold dust, atomic and molecular hydrogen sug- 
gests that when the gas is removed also the dust is pulled out 
the galactic disk. These results provide interesting insights into 
the evolution of an extreme case among perturbed cluster galax- 
ies. Once combined with the discov ery of truncated dust disks 
in Hi-deficient Virgo cluster spirals (ICortese et al.ll201()bl) . our 
analysis clearly highlights the great potential oi Herschel for our 
understanding of the effects of the cluster environment on the 
dust properties of galaxies. 
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